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Abstract 

A spherical shell model for turbulence, obtained by coupling N 
replicas of the Gledzer, Okhitani and Yamada shell model, is con¬ 
sidered. Conservation of energy and of an helicity-like invariant is 
imposed in the inviscid limit. In the A —> oo limit this model is 
analytically soluble and is remarkably similar to the random coupling 
model version of shell dynamics. We have studied numerically the con¬ 
vergence of the scaling exponents toward the value predicted by Kol¬ 
mogorov theory (K41). We have found that the rate of convergence to 
the K41 solution is linear in 1/N. The restoring of Kolmogorov law has 
been related to the behaviour of the probability distribution functions 
of the instantaneous scaling exponent. 
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The small-scale statistics of a fluid in a regime of fully-developed turbu¬ 
lence is an interesting and still open question. Among the many problems 
which have to be faced there is the understanding of anomalous scaling of 
structure functions Sp{r) defined as: 

Sp{r) =< \v{x -I- r) — v{x)\^ >=< \6rv\^ >~ r = |r|, (1) 

for r in the inertial range, i. e. for those scales, much smaller than the integral 
scale and much larger than the viscous scales, where inertial non-linear terms 
of Navier-Stokes (N-S) equation are dominant. 

The phenomenological theory of Kolmogorov (K41) is based on the hy¬ 
pothesis of an energy cascade from large to small spatial scales, local in scale 
size, in which all statistical information about large scale is lost, save for the 
mean energy-dissipation rate. It leads, by simple dimensional arguments, to 
the celebrated K41 law: (p^^ = p/3 |l|. This result is exact for p = 3 and in 
good agreement with the experimentally and numerically measured exponent 
for p = 2 - the energy spectrum scaling E{k) ~ is very well fitted - but 
there are strong evidences of its breakdown for larger p. 

Phenomenological theories that modify the Kolmogorov approach assign¬ 
ing a key role to the statistics of the spatial distribution of energy dissipation 
have been introduced in the last thirty years [0-^]. These approaches are 
able to fit with good accuracy the intermittent deviations from the p/3 law 
but do not have any direct link with the dynamical evolution of N-S equa¬ 
tions. On the other hand, the goal of analytically obtaining the corrections to 
K41 exponents directly from N-S equations is far from being achieved. Until 
now, neither in closure theories nor by a renormalization group approache 
have any satisfying results been obtained on this problem. 

A very interesting approach was proposed by Kraichnan in 1961 |^. In 
the so-called random coupling model (RCM) he replaced the N-S equation 
by a set of modified equations having the same basic structural properties 
of N-S: quadratic non-linearity, non-linear quadratic invariants, existence of 
truncated inviscid equipartition solutions. The velocity field is replaced by 
N fictitious random fields u°{k,t) satisfying the equations: 

Kik,t) = -^kmPijik) ( 2 ) 

where ^a/Sa are quenched random phases subject to suitable restrictions (see 
[^]). The closure problem can be solved in the N ^ oo limit and gives the 
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same equations of direct-interaction approximation (DIA), a two-point clo¬ 
sure theory introduced earlier by Kraichnan himself |^. Recently Mou and 
Welchman [Q introduced an alternative large-N model, the spherical model 
(SPM), by taking non-random <ha/3o- coefficients and imposing additional sym¬ 
metries on the system. Their idea is basically that of 1 /N expansion of critical 
phenomena: one generalizes the model under consideration to one with an 
higher symmetry (in standard application, for example, one considers 
the rotation group 0{N), with physical values N = 1,2,3); the N ^ oo 
limit is analytically soluble and a systematic expansion in powers of 1/A^ 
may be developed |^. Mou and Welchman, following the work of Amit and 
Roginsky [§, chose the group Gn = 0(3) for all N but allowed the dimension 
of the representation to diverge with N. The remarkable point is that the 
SPM and the RCM are similar: in both the inhnite N limit leads to the DIA 
equations [^. Both models, then, are plagued with the same problem: DIA 
equations violate the basic physical principle of galilean invariance and give, 
consequently, the wrong energy spectrum scaling, E{k) ~ xhe correct 

-5/3 exponent has been recovered, in fact, in a Lagrangian modihcation of 
DIA in which galilean invariance is restored (for an overview about these 
theories see 0). So it would be ill advised to try to get anomalous correc¬ 
tions of scaling exponents starting from the wrong zeroth order and taking 
successive orders in which the same problems occur. However, the basic idea 
of getting the intermittent corrections as asymptotic expansions in powers of 

1 /vn i-ini 


1 . e. 


5Cp-cr-c. = E 


k>l 


Sp) 

J^k/2 


(3) 


can be tested on simple models in which these inconsistencies do not show 
up like, for example, the shell models, as Eyink has recently proposed . 


Shell models are dynamical systems in which, exploiting the idea that 
active degrees of freedom in a turbulent flow are hierarchically organized, 
a reduced number of variables (one or two) are considered for each shell of 
wavevectors with modulus k G {kn,kn+i], where k^ = A:o2". Among shell 
models, the so-called GOY (from Gledzer, Ohkitani and Yamada) [^-[|T^ is 


of considerable recent interest. The model is dehned by the following ODE’s: 


d 

dt 


u. 


= -uk^Ur, 


+ ikr, 


'^n+l'^n+2 


(1 -£) 


Un-lUn-2 


+ /n<5i 


n,l 

(4) 
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with n = 1, ...,Nsheiis and boundary conditions with null m_i, uq, 
and UN^f^^ii^+ 2 - The coefficients of the non-linear terms have been chosen in 
order to keep the energy constant in the inviscid unforced case; besides, for 
£ = 1/2, one has an inviscid invariant H = |j^ that is a 

sort of “shell model helicity” similar to the helicity H = f(kx v{k)) ■v{k)dk 
conserved in the 3d Euler equations. Considering that should describe 
velocity fluctuations at scale r = the structure functions are dehned by: 




I Ur. 


> 


(5) 


and show, in the inertial range, a scaling behaviour: Sp ~ This model 
has the remarkable property of reproducing quantitatively the exponents (p 
experimentally measured when e = 1/2 

In [|I3 Eyink dehned and studied the inhnite N limit of the following 


complex spherical shell models (SSM): 



with Anmi local couplings, i.e. vanishing outside a hnite range of neighbors 
and given by (pil|). This kind of models reduce to real shell model 

when N = 1 and the DIA equations obtained in the inhnite N limit have 
stationary solutions with K41 scaling when a hxed input of energy by the 
external force is imposed. 

We considered a slightly modihed version that gives the original GOY 
model at Y = 1. We generalized the dynamical equations for the real and 
imaginary parts of the complex variable Un = Xn + iyn to the following 
dynamical system: 


[jt + 

[jt + "*^0 


+ l/n+Yn+2 “ ^{x^n-lVl+l + vl-lXl+i) 




{xi_iyl_2 + yn-lXl_2)\* + r^nl 




(7) 


- yn+iyl+2 “ ^{x^-lXl+i - 2/f_l2//+l) 


/3 ^7 


{x^-lXl_2 - yn-iyl-2)]* + r^nl 


( 8 ) 
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where sums over repeated indices are implied; a = — N = 2J+1 

and n = ...^Nshells- To get the original model in the limit J — 0 we 

have taken x" and ?/" complex variables satisfying = (—l)“(x")* and 
= (—1)"(2/")* and = (—!)“(/")*. So there are 2iV = 2(2J + 1) real 
degrees of freedom per each shell. Following we imposed the invariance 
of the system under the transformation: 


< - Dig{U)vl (9) 

where the are the Wigner matrices of quantum mechanics (N- 

dimensional irreducible representations of the group SU{2)) and n" = (x“, y") 
So one must take: 


wf’’ = (-1)“ < JJJ - a| JJ/37 > 


( 10 ) 


which are the Clebsh-Gordan coefficients. The properties of these coeffi¬ 
cients, together with the fact that we considered only even J’s, guarantee: 

(i) invariance of the system under any permutation of the replica indices; 

(ii) conservation of volume in phase-space; (iii) conservation, in the inviscid 
unforced case, of energy: 

1 N +J 

E Kt+Kr (11) 

^ n=l a=-J 

and of shell-model helicity: 

1 N +J 

^f = ;iE E (-i)”*:»ll<P + l»:tl' (12) 

^ n=l a=-J 

Let US stress at this point that the choice of considering the SSM rather 
then a shell RCM has been dictated by the considerable advantage that the 
former present in a numerical analysis, avoiding the averages on the random 
coupling coefficients. 

We analyzed the scaling properties of this model at varying the number 
of fields n" (up to a maximum number of 21 replicas). As the moments of 
the PDF’s of the velocity variables are affected by period three oscillations 
(responsible for large uncertainties in scaling exponents) and as quite a large 
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precision is needed to calculate deviations from K41 solution we considered 
the scaling properties of the following suitable triple products: 


Un{N) = -n] A„+i(7V) + 




(13) 


where: 


A„(A) = 
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ap'y 

+2/n (• 


X, 


f P 7 I P 7 
(<+l2/n-l + Vn+lK-l 


P 


n+l^n—1 yn+lVn—l 


(14) 


The variables Iln, have a physical relevance, Unkn representing the energy 
flux from the scale n to larger wavenumber scales, and are not plagued with 
periodic oscillations (see £g. 1). We make use of a relation analogous to the 
structure equation of Kolomogorov as in |]T^ and so the following scaling 
law is expected: 

S”(/V) = (|n„(Af)y3)~l;„-«''> (15) 

In hgure 2 we plotted the behaviour of 6(p{N) = — Cp(^) as function 

of 1/N for different values of p. The corrections to the K41 solution tend to 
zero when A —>■ cx) and the rate of convergence is linear in 1/N. It should 
then be possible to develop an asymptotic expansion for the ^Cp’s of the form 


as Eyink suggested for the shell models in |10 


In our numerical simulations we used Nsh^us = 15 and 19, with u = 10“® 
and 10“®, /" = 5 10“^(1 +i), fco = 6.25 10“^ and a number of replicas up to 
21. Equations have been integrated by a slaved Adams-Bashforth algorithm 
and temporal averages are over several thousands of eddy turn-over times of 
the hrst shell. 

Beyond these numerical results an interesting physical problem is to un¬ 
derstand the dynamical mechanism by which K41 is restored. As Kraich- 
nan pointed out in |^, the K41 law is not ruled out just because energy- 
dissipation fluctuates: it could be a priori possible if there were spacewise 
diffusion effects of sufficient strength to suppress fluctuations of energy trans¬ 
fer in the inertial range. We suppose that a similar mechanism is present 
in the SSM, although it is not completely clear to us to what extent the 
similarity can be pushed. In the GOY model the anomalous scaling law is 
accompanied by an intermittent behaviour of the energy transfer in the in¬ 
ertial range with very large fluctuations. A quantity that measure the local 
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(in time) singularity is the instantaneous scaling exponent for fluxes, defined 
as: 


hit) 




+ 1 


£ 1^2 


n. 


n, 


n+l 


( 16 ) 


where rimin and Hmax delimit the inertial range. In a laminar regime one has 
h = 3, h = 1 corresponds to a K41 scaling while a value of h < 1 represents 
a more singular velocity held. Such a value of h is realized during the fast 
energy bursts [0. We measured h(t) and calculated its PDF for different 
values of N. We observed that increasing the number of coupled replicas the 
peak of the distribution gets sharper and the minimum value of h increases, 
indicating that the held is getting less singular in the inertial range (see hg. 
3). This result could be due to possible exchanges of energy between the 
increasing number of degrees of freedom at the same scale. 

To conclude we would like to stress that the SSM is the only model 
we know in which it has been shown that there is approach to the K41 
solution by tuning a parameter (the number of replicas) starting from an 
intermittent solution. Let us remark, anyway, that in a diherent contest, 
that of passive scalar advection, two groups have been succeded in calculate 
anomalous scaling exponents by expansion around limiting gaussian cases 
[O, M. At now nothing similar has been done with the N-S equation also 


if it is believed that an inhnite-dimension limit should converge to K41. 
Beyond the possibility of getting an asymptotic expansion for the anomalous 
corrections of the form (^) we believe that further work in this direction 
should deepen the comprehension of intermittency in real turbulence. 
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FIGURE CAPTIONS 


• FIGURE 1: 

Structure function Sp^g(A^) vs kn in a log-log plot for = 5 and 
Nshells = 19 (crosses with line). We plotted also the scaling behaviour 
of the 6-th moment for a = 1 (diamonds). 

• FIGURE 2: 

Anomalous corrections 6(p{N) versus 1/N for p = 5,..., 11- The error 
bars have been obtained by a least square £t. Note that the origin is 
an analytically calculated point. 

• FIGURE 3: 

Probability distribution function of the instantaneous scaling exponents 
h(t) of equation (jT^)for different values of N. 
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